NASA/TP-1998-208707 



Track Structure Model for Radial 
Distributions of Electron Spectra and 
Event Spectra From High-Energy Ions 

F. A. Cucinotta 

Langley Research Center, Hampton, Virginia 
R. Katz 

University of Nebraska, Lincoln, Nebraska 
J. W. Wilson 

Langley Research Center, Hampton, Virginia 


National Aeronautics and 
Space Administration 

Langley Research Center 
Hampton, Virginia 23681-2199 


September 1998 



Available from the following: 


NASA Center for AeroSpace Information (CASI) 
7121 Standard Drive 
Hanover, MD 21076-1320 
(301) 621-0390 


National Technical Information Service (NTIS) 
5285 Port Royal Road 
Springfield, VA 22161-2171 
(703) 487-4650 


Abstract 


An analytic method is described for evaluating the average radial electron 
spectrum and the radial and total frequency-event spectrum for high-energy ions. 
For high-energy ions , indirect events make important contributions to frequency- 
event spectra . The method used for evaluating indirect events is to fold the radial 
electron spectrum with measured frequency-event spectrum for photons or elec- 
trons. The contribution from direct events is treated using a spatially restricted 
linear energy transfer (LET). We find that high-energy heavy ions have a signifi- 
cantly reduced frequency-averaged lineal energy (y F J compared to LET \ while 
relativistic protons have a significantly increased y F and dose-averaged lineal 
energy (y for typical site sizes used in tissue equivalent proportional counters. 
Such differences represent important factors in evaluating event spectra with lab- 
oratory beams , in spaceflight , or in atmospheric radiation studies and in valida- 
tion of radiation transport codes. The inadequacy of LET as descriptor because of 
deviations in values of physical quantities , such as track width , secondary elec- 
tron spectrum and y^for ions of identical LET is also discussed. 


Introduction 

This paper presents an analytic model for describ- 
ing the radial distribution of electrons and frequency- 
event spectra measured by proportional counters used 
in the dosimetry of radiation fields in space and the 
upper atmosphere (refs. 1 and 2). At high energies the 
track width of an ion will extend to 100’s of microns 
or more because of delta-ray diffusion. Monte-Carlo 
models have difficulty with these calculations because 
of large computational times and poor statistics in con- 
sidering events at distances greater than 1 pm from the 
track of an ion (ref. 3). Furthermore, in most applica- 
tions with high-energy nuclei including space or atmo- 
spheric radiation studies, nuclear fragmentation and 
energy loss processes will lead to a broad spectrum of 
ion types and velocities. These considerations point to 
the usefulness of analytic approaches to treat track 
structure for high-energy particles. An analytic model 
is developed to evaluate frequency-event spectra that 
include the indirect events important for high-energy 
ions. Furthermore, the radial dependence of the spec- 
trum is described. The model is being used to compare 
the results of radiation transport codes directly to 
recent spaceflight measurements as described else- 
where. These comparisons will provide important val- 
idation of radiation transport codes and environmental 
models. 

Track structure models describe the relationship 
between the spatial distribution of energy deposition 
in the form of positions of ionization and excitation. 


the geometric structure of target molecules, and the 
response of a physical or biological system. The first 
track model formulated was the average-track model, 
or amorphous-track model, which considered the 
radial dose about the path of the ion as the descriptive 
parameter for ion effects. This model has been used 
for over 30 years in describing the response of physi- 
cal and biological systems (refs. 4-6). The average- 
track model has the advantage of simplicity; however, 
the model considers only simplified target geometry 
and ignores fluctuations in energy deposition. Monte- 
Carlo track simulations consider local fluctuations in 
energy deposition and provide methods for dealing 
with complicated target geometry including applica- 
tions to treat details of DNA molecular structure and 
DNA folding (refs. 7 and 8). However, the average- 
track model remains the most successful parametric 
approach for describing the response of biological sys- 
tems and physical detectors to ions. 

The importance of track width has been discussed 
for many years in relationship to radiation quality and 
the resulting inadequate capability of linear energy 
transfer (LET) or the ratio of the square of the effec- 
tive charge to ion velocity (Z* 2 /|3 2 ) as indicators of 
biological effectiveness (ref. 9). The parameter Z* 2 /| 3 2 
was noted by Katz to be useful only when comparing 
ions of similar velocity (ref. 9). Experimental observa- 
tion of the importance of track width has been demon- 
strated for endpoints such as inactivation of yeast, bac- 
terial and V79 mammalian cells, and mutation of the 
HPRT gene in V79 cells (refs. 10-12). Recently, 



calculations with the average-track model (ref. 6) 
showed that, for cellular mutation, the spatial distribu- 
tion of sites for cellular inactivation and gene mutation 
necessitates the use of a track structure description for 
radiation quality that indicates the inadequacy of the 
parameters LET or Z* 2 /p 2 . Our calculation of the 
radial distribution further illustrates inadequacies in 
these parameters when biological effectiveness is 
dependent on electron energy. 

Algorithms for converting particle energy spectra 
to lineal energy spectra are needed to correlate radia- 
tion transport code predictions to lineal energy spectra 
measured by proportional counters. Microdosimetric 
approaches to energy deposition have relied mostly on 
Monte-Carlo simulations (refs. 13 and 14), which are 
computationally inadequate for describing diverse 
radiation fields such as that seen in space or in the 
upper atmosphere. Several factors are known to be 
important in evaluating lineal energy spectra. For 
high-energy ions, a significant fraction of events occur 
from secondary electrons produced by ions that do not 
directly enter the volume of the detector. Other factors 
include the treatment of secondary electrons trans- 
ported outside the volume by ions that do pass through 
the volume, nuclear reaction effects, and the effects of 
straggling of the ions in the volume. For each of these 
factors the role of wall composition and thickness 
must be studied. Also, for low-energy ions there is a 
significant change in the rate of energy loss inside the 
volume, including the effects of stopping ions 
(ref. 15). 

This report presents an approach for treating the 
indirect events from delta rays and introduces a spa- 
tially restricted energy deposition model for direct 
events. The other factors noted will be discussed 
elsewhere. 

Model for Radial Electron Spectrum 

The approach of the average-track model has been 
to consider the primary electron spectrum from ion 
interactions with target atoms and fold this spectrum 
with average transmission properties of electrons to 
obtain the spatial distribution of electron dose as a 
function of radial distance from the path of the ion. As 
introduced by Kobetich and Katz (ref. 4) the radial 
dose is given by 


D « ( ' } = 2^ X J dn J d(0 J7 [E <*> “W'* “)] 

Xn(,,a>) (dCD (1) 

In equation (1) £ is the residual energy of an electron 
co after travelling distance t, and q(f,co) is the transmis- 
sion probability that an electron with starting energy co 
penetrates a depth t. We have included an angular dis- 
tribution for the primary electrons with energy co and 
solid angle ft. The subscript 5 indicates that it is the 
dose contribution from ionization by secondary elec- 
trons at a radial distance t from the path of the ion. The 
input functions for the evaluation of equation (1) are 
described by Cucinotta et al. (ref. 16). Note that the 
cross sections for electron production from protons are 
scaled to heavy ions using effective charge. 

The LET can be described by integrating the 
radial dose distribution over all radial distances and 
including other contributions such as excitations, 
nuclear stopping, and high-energy corrections (ref. 1 7) 

LET = 2ti f M t dt[D & (t ) + D exc (f)] 

+ Nuclear stopping 
+ High-energy corrections (2) 

where / ^ is the maximum distance of electron penetra- 
tion. Brandt and Ritchie (ref. 18) have considered a 
formulation of the excitation term, D exc (f), as 

D excW = C exc exp(-//2d)/f 2 (3) 

with d = p/2 tOj. with co r = 13 eV for water. The radially 
restricted LET can be introduced by limiting the upper 
limit in equation (2) as defined by 

L r = 2 ti ft dt[D s (t) + D exc (t)] (4) 

Equation (4) ignores the nuclear stopping and 
high-energy corrections to the LET. In many applica- 
tions the number of electrons, as well as their energy 
spectrum, is required for describing the response of a 
system. The average or residual energy spectrum of 
electrons penetrating to a radial distance t is derived 
from equation (1) as 
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( 5 ) 


where the primary electron energy oo is now a function 
of the residual energy. Equation (5) shows that the 
attenuation of the electron spectrum is through two 
factors: first an overall factor of Mt for all secondary 
electrons, and second an additional attenuation for 
low-energy electrons dependent on their starting 
energy and depth of penetration. 

Figures 1 and 2 show calculations of the radial 
dose for ions of linear energy transfer (LET) of 30 and 
150 keV/pm. These results show the large differences 
in radial energy deposited due to differences in track 
width, which is dependent on ion velocity. Such dif- 
ferences are expected to be important for specific 
target molecule sizes and the spatial distribution of 
these molecules. Also shown in figure 1 are results 
based on an assumption that the electrons are ejected 
normal to the ion path leading to an underestimate of 
the dose at small distances. Figures 3 and 4 show the 
secondary electron spectrum for these same ions. We 
have plotted the spectrum from equation (5) versus 
energy for several impact parameters and also we plot 
2 ntE <[>(£, t) versus energy on a linear-log scale to 
show the fractional contribution from each decade of 
electron energy. Lower charge ions are seen to have 
both a confined track width and electrons of lower 
energy (i.e., electrons of higher LET). Experiments 
with soft X rays indicate increased biological effec- 
tiveness for electrons with energies less than several 
keV (ref. 19). 

For ions of a given value of LET, low charge and 
energy ions (LZE) may have increased effectiveness 
because the predominance of low-energy electrons in 
comparison to high charge and energy ions (HZE). 
However, for large target volumes, including a 
response dependent on alterations in spatially distrib- 
uted target molecules, the present energy deposition 
model would predict that the effects of track width 
would increase the effectiveness of HZE ions in com- 
parison to LZE ions. The secondary electron spectrum 
varies more slowly as a function of radial distance for 
the LZE ions and contains a larger fraction of low- 
energy electrons both at small radial distances and at 


the maximal radial distances where the highest energy 
ejected electrons are stopped (electron track ends). 
Figures 3 and 4 show that for ions with atomic num- 
bers less than iron (Z < 26), overlapping electron 
tracks in small volumes such as the nucleosome would 
occur only at very close distances to the path of the ion 
(<10 nm). 

Direct Events 

For fast charged ions undergoing small velocity 
changes in the site, the energy deposited for path 
length x is given in terms of the linear energy transfer 
as 

e = Lx (6) 

For high-energy ions a significant fraction of the 
LET escapes the volume of a site by electron transport 
and therefore does not contribute to the events inside 
the site. The energy deposited, as restricted by the spa- 
tial distribution of secondary electrons, is given by 

e = J dx'L r (x') ^ 

For a given path length the radial distance from 
the ion to the center of the volume is 



where d is the diameter of the site. For a spherical site 
the relation given by equation (8) limits the radial 
extension of the track for direct events such that the 
maximum radial extension of the track at a polar angle 
9 is given by 

t M (x, (?) = 2/y t x cos <p (9) 

where 

r x , = x' for x' < d/2 (10a) 

r x , = d-x' for x' > dll (10b) 

The energy deposited for a given path length as 
restricted by electron diffusion is then given by 

r* r <P) 

£ = )dx'jd(?j tdtD(t) (11) 

For a distribution of path lengths^), the distribu- 
tion of event sizes is given by 
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/(e) dz = Nf(x) dx 


( 12 ) 

Such that 

/(e) = A/[e(oO](^J (13) 

where from equation (11) 

— = 1^1 tdt D(t ) (14) 

The lineal energy is given by y = z / x\ where x A is the 
average path length in the volume and the distribution 
in lineal energy is given by 

f(y ) = Nx A f[y(x)]f^y (15) 

Low velocity ions with insufficient range to transverse 
the site at a given path length will deposit all their 
energy in the site and the direct events should be 
appropriately corrected (ref. 15). The distributions for 
direct events should also be corrected for ion strag- 
gling and nuclear reaction effects should be described. 

Indirect Events 

Particles that do not pass directly through the vol- 
ume deposit energy through secondary electrons; these 
events are denoted as indirect events. Figure 5 shows 
the fraction of the LET from iron ions as a function of 
their energy that are from indirect events in a homoge- 
neous medium for several site sizes. At high energies, 
a significant fraction of the LET from indirect events 
occurs for site sizes in the 0.5-4 Jim range, which is 
used frequently for radiation field definition. For small 
site sizes expected to be important in producing DNA 
damage and mutations, the LET is dominated by the 
so-called indirect events. 

The average distribution of secondary electrons as 
a function of radial distance can be used to evaluate 
the contribution from indirect events using a folding 
approach. We treat this spectrum as an isotropic 
source incident on the sensitive volume. The theoreti- 
cal evaluation of microdosimetric spectra from elec- 
trons is difficult to treat analytically because of the 
small mean free path of electrons for both elastic and 
inelastic collisions and the importance of energy and 
range straggling for electrons. Extensive measure- 
ments of event spectra using photons and in some 
cases electrons over a large range of energies have 


been made with TEPC’s (refs. 20—22). For smaller site 
sizes, Monte-Carlo calculations have been made for 
electrons with energies from 0.1 to 100 keV (refs. 6 
and 7). Our approach is to consider the average elec- 
tron spectrum of electrons at radial distance t from the 
path of the ion and fold this distribution with represen- 
tations of measurements for photons or electrons. In 
this way, local fluctuations in energy deposition from 
delta rays are taken into account. 

The event spectrum from indirect events is found 
by assuming the radial distribution of electrons is inci- 
dent on a spherical site at distance t from the track cen- 
ter and folding this spectrum with the event spectrum 
for electrons of a given energy 

•^out^’ = J dE <j)(t, £)/ e ( E,y ) ( 16 ) 

where the distributions / e _ (£, y) are the events by elec- 
trons of energy £ in a particular volume as inferred 
from experiments (refs. 20-22). 

The y spectra from measurements with photons 
are parameterized as 

f(E,y) = w|cexp^j+ (1 -c) exp (17) 

with N as a normalization constant and the parameters 
chosen as b = 6.5 keV/|am, and 


a = + a 2 exp 


E \] 

photon \ 

60 J 


(18) 


<«> 

Values for the parameters a h a 2 , and c 0 are given in 
table 1 for several spherical sites with diameters from 
0.5— 4.0 fim. Equation (17) provides a good approxi- 
mation to measured values for vp and y D as a function 
of photon energy as shown in figure 6. For an expo- 
nential spectrum, the relationship yp/yp = 2 is found 
and is observed in most experiments below 100 keV. 
The photon spectrum of equation (17) is fit to 
experiments for photon irradiation in walled counters. 
Some differences in the response of walled and wall- 
less counters for photons and electrons should be 
expected and are not described here. To relate the 
photon energy to its secondary electrons we use the 
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average secondary electron energy from photon 
irradiation resulting from Compton scattering and the 
photoelectric effect. Event spectrum with photons for 
wall-less counters indicates about a 10- to 20-percent 
reduction in y F in comparison to walled counters. 
Also, measurements made directly with electron 
beams (refs. 21 and 22) suggest slightly lower values 
for y F and y D than those derived here from photon 
exposures. The modification of equation (18) to 
^wall-less * 0.8a wall will approximately account for 
some of these differences; however, such effects need 
to be studied further. The total event spectrum is found 
by integrating overall radial distances and including 
the contribution from direct events. 


Table 1 . Parameters for Photon Lineal Energy Distribution 


Site diameter, 
pm 


a 2 

^0 

0.5 

0.3 

4.6 

0.0050 

1.0 

0.14 

4.2 

0.0015 

2.0 

0.14 

3.4 

0.0010 

4.0 

0.18 

2.6 

0.0005 


Results and Discussion 

Figure 7(a) shows the frequency average of the 
specific energy, zp(t) as a function of radial distance 
from the path of 600 MeV/amu iron ion for a site of 
1.3 pm diameter. The lineal energy spectrum 
described above is converted to specific energy using 
the relationship 

Z = 0.204 ^ (20) 

d 

where z is in Gy, y is in keV/pm, and d is the site 
diameter in (am. Good agreement with the experiment 
of Metting et al. (ref. 23) is found. The frequency- 
averaged values and their correlation with known val- 
ues for electrons and photons clearly indicate the role 
of delta rays in outside events. Figure 7(b) shows com- 
parisons of calculation to experiment for the mean 
square of the specific energy, which is given by 

z 2 D (t) = z D (t) D(t) + D 2 (t) (21) 

Again good agreement between experiment and model 
is found. Figure 8 shows comparisons of model to 
experiments (ref. 24) for z D (t) for 13.7 MeV/amu ger- 
manium ions in 0.5 and 1 .0 pm sites. Comparison of 


figures 7 and 8 indicates a harder electron spectrum 
for lower energy ions. The comparisons of the radial 
distribution for average specific energy provide new 
support for the success of the average-track model of 
Katz (refs. 5 and 9) in describing relative biological 
effectiveness for diverse radiation fields. These com- 
parisons indicate that models that describe local fluc- 
tuations in energy deposition would account for 
factors responsible for the small differences in 
response seen between gamma-rays and hard X rays. 
For a parametric model summing the effects of expan- 
sive electron spectra such differences may be of 
reduced importance; however, such differences could 
be described in the present approach. 

Figures 9-1 1 show calculations of the frequency- 
and dose-averaged linear energy versus kinetic energy 
for proton, oxygen, and iron beams. Also shown are 
values of the LET. High-energy proton beams are seen 
to have much higher values for y F and y D than LET, 
with y D exceeding LET by a factor of about 5 at high 
energies in a 1 pm site. Such difference will have large 
impact in validating radiation transport codes (refs. 25 
and 26) with TEPC measurement (ref. 2). For heavy 
ions, values of y^y agree with LET within 20 percent; 
however, these ions carry an additional low y compo- 
nent to their spectrum as shown in figure 12. Such 
observations were seen in the experiments of Dicello, 
Wasiolek, and Zaider (ref. 27) and the present model 
provides an analytic model to describe these effects. 
For studying energy deposition from relativistic ions, 
the use of large site sizes (>1 pm) would not signifi- 
cantly diminish the contribution of indirect events. For 
space radiation studies, large site sizes would further 
distort contributions from low velocity ions produced 
through target fragmentation (ref. 26). 

Conclusions 

The average-track model has been extended to 
describe the radial distribution of electrons about the 
path of an ion. We have used this spectrum in a fold- 
ing model to predict the contribution of indirect events 
to frequency-event spectra. In this method, local fluc- 
tuations in energy deposition are included in the 
model. Tissue equivalent proportional counters are 
often used to measure the spectrum of radiation types 
onboard spacecraft or in the upper atmosphere. The 
method used here will provide response functions to 
validate the results of radiation transport codes and 
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environmental models with such measurements. Large 
differences between linear energy and lineal energy 
spectra for relativistic ions are expected, based on the 
model described, because of indirect events. Most 
importantly, high-energy light ions have significantly 
increased values for frequency- and dose-averaged lin- 
eal energies compared to LET. Finally, we have dis- 
cussed differences in track width, secondary electron 
spectrum, and frequency event spectrum for ions of 
identical LET. Clearly these differences, along with 
other noted track structure effects, preclude the use of 
LET or Z* 2 / p 2 as descriptors of radiation quality. 

NASA Langley Research Center 
Hampton, VA 23681-2199 
April 6, 1998 
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Dose, Gy 




(a) For 1 MeV protons. (b) For 377 MeV/amu neon. 

Figure 1 . Comparison of calculations to experiments for radial dose for ions of LET = 30 KeV/jim. 




(a) For 0.75 MeV/amu alpha particles. (b) For 1000 MeV/amu iron particles. 

Figure 2. Calculations of radial dose distributions from ions of LET of 1 50 KeV/|xm. 
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Electron energy, E, keV Electron energy, E, keV 


(a) Electron spectra for 1 MeV protons. (b) Electron lethargy spectra for 1 MeV protons. 



Electron energy, E, keV Electron energy, E, keV 


(c) Electron spectra for 377 MeV/amu neon. (d) Electron lethayrgy spectra for 377 MeV neon. 

Figure 3. Calculations of electron energy ( E) spectrum <p(£, t) or 2 nt <p (£, t) at various distances / from ions of LET of 
30 keV/pm. 
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(c) Electron spectra for 1000 MeV/amu iron ions. (d) Electron lethargy spectra for 1000 MeV/amu iron ions. 

Figure 4. Calculations of electron energy (£) spectrum q> (£, r) or 2nt <p(£, /) at various distances t from ions of I F T of 
150 keV/pm. 
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Photon energy, kV 
(c) Diameter = 2 fim. 



yp - fit (diameter = 4.0 |im) 

>D - f,t 

yp - experiment (diameter = 4.0 pm) 
yp experiment 





D (t) and z D z (t), Gy Zp (t), Gy 



(a) Average specific energy, zp(0- 



(b) Square of mean specific energy. 


Figure 7. Comparison of calculation to experiments (ref. 23) for average energy in 1.3 pm site as function of radial distance t 
for 600 MeV/amu iron ions. 
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Distance, t, Jim 


Figure 8. Comparison of calculation to experiment (ref. 24) for average specific energy in 0.5 and 1.0 pm sites as function of 
radial distance for 13.7 MeV/amu germanium ions. 
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or LET, keV/jim 



Figure 10. Calculations of dose-averaged lineal energy (y D ) versus kinetic energy for oxygen ions in 1 |am sites. 
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(a) Helium. 



(b) Iron. 

Figure 12. Comparison of event distributions, yf(y), in 1 pm site for helium and iron ions of 200 MeV/amu. 
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